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Abstract 

We report first observations of the decays B~ — > D s K~tt~ , using 292 fb~ of data collected at 
the T(4S) resonance energy by the BABAR detector at the PEP-II e + e~ collider. The branching 
fractions are measured to be B{B~ -> DfK~ir~) = (1.88 ± 0.13 ± 0.41) • 10~ 4 and B{B- -> 
D* + K-tt~) = (1.84 ± 0.19 ± 0.40) • 10" 4 . 
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1 INTRODUCTION 



First evidence for so-called inclusive flavor correlated production of Df in B~ decays was reported 
recently jTj with a branching fraction of B(B~ —* DfX) = (1.2 ± 0.4)% These decays are 
mediated by a 6 — > c quark transition and require at least three final state particles, including the 
production of an ss pair from the vacuum ( ss "popping"). An example for a three-body B~ decay 
with a Df in the final state is B~ — > Df K~ir~ . The corresponding B° decay is B° — > DfK°7t~. 
The Feynman diagram for B~ — > Dg*^ + K~tt~ decays is shown in Fig.^ In case of B° — > DfK°7t~ , 
an additonal contribution from a VF-exchange diagram with ss and (id popping may exist. If we 



d n 




Figure 1: Feynman diagram for B — > DfK ir . 

replace the 7r _ in Fig. ^ which comes from the hadronization of the W~ boson with a K~ , we 
get the Cabibbo-suppressed decays B~ — > DfK~K~ and B° — > Df K° K~ . It is interesting to 
note that the final state DfK°K~ can also be reached from a B° decay. In this case the decay is 
mediated by a b — > u quark transition, but the hadronization is not Cabibbo-suppressed. Thus 
a B° can either decay directly to DfK°K~ or via mixing followed by B° — > DfK°K . The 

interference between the two decay amplitudes for decay with and without B°B° mixing leads to 
a time-dependent CP-asymmetry that is sensitive to sin(2/J + 7). In case the contribution from 
the higher D** resonances decaying into DfK turns out to be large, it may also be interesting to 
measure the resonant parameters independently from the analysis using B — > Dtttt decays [Sj. 
No exclusive B~ — > D s X or B u — > Dg X decay mode has hitherto been observed. Limits 

on the branching fractions from the analyses by other experiments are listed in Table Q In this 

— (*)+ ~ — 

paper we report the first measurement of the decay modes B — > D s K it . 

2 THE BABAR DETECTOR AND DATASET 

The analysis uses a sample of approximately 292 fb _1 , which corresponds to about 324 million 
T(4S) decays into BB pairs collected with the BABAR detector at the PEP-II |lj asymmetric- 
energy B-factory. The BABAR detector is described elsewhere and only the components crucial 
to this analysis are summarized here. Charged particle tracking is provided by a five-layer silicon 
vertex tracker (SVT) and a 40-layer drift chamber (DCH). For charged-particle identification, 
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Table 1: Upper limits from ARGUS i(>J and CLEO [Tj on B~ -> L>r )+ ET"7r~ branching fractions. 



Experiment 


Decay Mode 


Upper limit (@90% C.L.) 


ARGUS 


B~ - 


-> D+K--K- 


8 x 10~ 4 




B- - 


-» DI+K--K- 


12 x 10~ 4 


CLEO 


B- - 


-> D+K--K- 


5 x 10~ 4 




B~ - 


-> D* + K--k- 


6.8 x 1(T 4 



ionization energy loss {dE/dx) in the DCH and SVT, and Cherenkov radiation detected in a ring- 
imaging device are used. Photons are identified and measured using a thallium-doped Csl-crystal 
electromagnetic calorimeter. These systems are located inside a 1.5 T solenoidal superconducting 
magnet. We use GEANT4 [S] software to simulate interactions of particles traversing the BABAR 
detector, taking into account the varying detector conditions and beam backgrounds. 

3 ANALYSIS METHOD 

The optimal selection criteria as well as the probability density distributions of selection variables 
are determined by a blind analysis based on Monte Carlo (MC) simulation of both signal and 
background. For the calculation of the expected signal yield we assume B(B~ -> D { s * )+ K-ir-) to 
be 10 -4 (i.e. about 10% of the measured B(B~ — > D + ir~ir~) [Hj). We use MC samples of our signal 
modes and, to simulate background, inclusive samples of B + B~ (784 fb _1 ), B°B° (774 fb _1 ), cc 
(247 fb _1 ), and qq, q = u,d,s (246 fb _1 ). In addition, we use large samples of simulated events of 
rare background modes which have final states similar to the signal. We have verified that our MC 
correctly describes the data by comparing distributions of various selection variables. 

Candidates for Df mesons are reconstructed in the modes — ► 4>tt + , K*°K + , and K®K + , 
with 4> K + K~ , K*° — » K~ir + and K® — ► 7r + 7r~ . The K® candidates are reconstructed from two 
oppositely-charged tracks, that come from a common vertex displaced from the e + e~ interaction 
point. We require the significance of this displacement (measured flight distance divided by an 
estimated error) to exceed 2. All other tracks are required to originate less than 1.5 cm away from 
the e + e~ interaction point in the transverse plane and less than 10 cm along the beam axis. Charged 
kaon candidates must satisfy kaon identification criteria that are typically around 92% efficient, 
depending on momentum and polar angle, and have a pion misidentification rate at the 5% level. 
The <j> — > K + K~ , K*° — > K~tt + and — > tt + it~ candidates are required to have invariant masses 
close to their nominal masses (we require the absolute differences between their measured masses 
and the nominal values [3] to be in the range ±15 MeV, ±50 MeV and ±10 MeV, respectively). 
The polarizations of the K*° and <fi mesons in the Df decays are employed to reject backgrounds 
through the use of the helicity angle 6h, defined as the angle between the K~ momentum vector 
and the direction of flight of the Df in the K*° or eft rest frame. The K*° and <fi candidates are 
required to have | cos 9h\ greater than 0.5. 

The D* + candidates are reconstructed in the mode D* + — > 7. The photons are accepted if 
their energy is greater than 100 MeV. The Df and D* + candidates are required to have invariant 
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masses in the interval [-10, 10] MeV/c 2 (for D+) and [-15, 10] MeV/c 2 (for D*+) from their nominal 
values [S] (the Df mass resolution is around 6 MeV/c 2 , and the asymmetric mass cut on D* + has 
an efficiency of about 90%). All Dg candidates are mass-constrained. The invariant mass of the 
D* + is calculated after a mass constraint on the daughter Df has been applied. Subsequently, all 
D* + candidates are subjected to a mass-constrained fit. 

We also require that photons from D* + are inconsistent with 7r° hypothesis when combined 
with any other photon having an energy greater than 150 MeV in the event (the 7r° veto window is 
±10 MeV/c 2 ). Finally, the B~ meson candidates are formed using the reconstructed combinations 
of D+K-tt- and D*+K~tt~. 

The background from continuum qq production (where q = u,d, s, c) is suppressed based on the 
event topology. We calculate the angle (6t) between the thrust axis of the B meson candidate and 
the thrust axis of all other particles in the event in the center-of-mass frame (cm.). In this frame, 
BB pairs are produced approximately at rest and have a uniform cos 6t distribution. In contrast, 
qq pairs are produced in the cm. frame with high momenta, which results in a | cos 9t\ distribution 
peaking at 1. | cos#t| is required to be smaller than 0.8. In addition, the ratio of the second and 
zeroth order Fox- Wolfram moments [10] must be less than 0.3. 

We extract the signal using the kinematical variables tties = \J E^ 2 — P*) 2 and AE = 

Si \J 771/2 + Pi 2 — wnere 1S the beam energy in the cm. frame, p* is the cm. momentum of 
the daughter particle i of the B~ meson candidate, and rrii is the mass hypothesis for particle i. 
For signal events, tiies peaks at the B~ meson mass with a resolution of about 2.6 MeV/c 2 and 
AE peaks near zero with a resolution of 13 MeV, indicating that the B~ candidate has a total 
energy consistent with the beam energy in the cm. frame. The B~ candidates are required to have 
\AE\ < 25 MeV (around 2a of the signal AE resolution) and m-Es > 5.2 GeV/c 2 . 

The fraction of events with multiple B~ candidates is estimated using the MC simulation and 
found to be around 3% for DfK~TT~ and 9% for D* + K~ir~ combinations. In each event with 
more than one B~ candidate that passes the selection requirements, we select the one with the 
lowest \AE\ value. 

After all selection criteria are applied, we estimate the B~ reconstruction efficiencies, excluding 
the subsequent branching fractions (see Table |2j). 

Table 2: Reconstruction efficiencies for B~ — > K~tt~ decays (excluding the subsequent 

branching fractions). 

Decay mode D+ <pir+ Df -> K*°K+ D+ -> K°K + 

B- -» DfK-n- 11.0% 7.0% 10.0% 

B- -» D*+K-7C~ 5.3% 3.4% 4.8% 



Background events that pass these selection criteria are represented by approximately equal 
amounts of qq continuum and BB events. We parametrize their ttt-es distribution by a threshold 
function [TT] : 

/(m E s) ~ rriEsVl ~ x 2 exp[-£(l - x 2 )], 
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where x = 2mEs/y / s, y/s is the total energy of the beams in their center of mass frame, and £ is 
the fit parameter. 

A study using simulated B and B + decays shows that some background events with distri- 
butions in tbes and in AE peaking near the signal region are expected in reconstructed B~ — > 
DfK~n~ candidates due to charmless and charmonium B~ decays with the same set of particles in 
the final state. For B~ — > D* + K~tt~ , no background of this kind is expected, due to the presence of 
the 7, which suppresses charmless and charmonium decay contributions. The peaking contribution 
is evaluated using the data by reconstructing ll D s " K~tt~~ combinations, where "Dj" candidates 
are selected from [±40, ±25] MeV sidebands around the nominal mass. In this procedure, we 
use the same selection requirements, as for the signal, except that candidates are not mass 

constrained. The resulting mES spectra are shown in Figure |^1 We fit the distributions using an 
extended unbinned maximum likelihood (ML) fit with a sum of a Gaussian (with a width and 
central value fixed from the MC simulation) and a threshold function /(otes) with the floating 
shape and normalization (see detailed expression of the likelihood function is Section EJ) . The fit 
yields 34 ± 12 events in the "signal" m ES peak for u D+"K-tt~ and 3 ± 7 for U D* + " K~~ir~ . Since 
the sideband interval is 1.5 times larger than the mass region used for signal selection, this 
translates into 23 ± 8 peaking background events expected in B~ — > DfK~ir~ . 
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Figure 2: tues spectra for the data "Df^K^ir^ (left) and K~n~ (right) combinations with 

no mass constraint applied on candidates, using Df mass sidebands [±25, ±40] MeV (1.5 

times the signal interval). 

We also study cross-feed between the signal modes and other decays with final states similar 

(*)+ — — — 
to our signal modes, including D s K~K~. The cut on AE of the B~ candidates effectively 

suppresses the cross-feed contributions, which do not exceed 2% of the reconstructed signal after 

all the selection criteria are applied. 

4 SYSTEMATIC STUDIES 

The summary of the systematic uncertainties is presented in Table |21 The total relative systematic 
error is estimated to be 22% for each B~ decay mode, with the largest contribution coming from 
the .D+ branching fractions uncertainty (15%). Other significant sources of systematic errors are 
found to be due to the difference between the selection efficiency in MC events and in the data 
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(estimated using the control mode B —* D s D°, D° — ► i^T 7r + ), and also due to the efficiency 

(*)+ — 

dependence on the D. s K invariant mass and its potential effect if the resonant contribution is 
present. 

Table 3: Summary of relative systematic errors (in %) for B- -» D [ * )+ K--k~ decays. 



Source B 


-* DfK 7r B 


-> D* + K 7T 


B counting 


1.1 


1.1 


MC statistics 


0.8 


1.4 


Tracking 


5 


5 


Particle identification efficiency 


4 


4 


-K^ efficiency 


0.5 


0.5 


7 (from D* s + — > Dg"f) efficiency 




2 


£> of sub-decays 


15 


15 


Peaking background contribution 


6 


3 


Cross-feed contribution 


1 


2 


Selection efficiency, Data/MC 


12 


12 


Signal and background shape uncertainty 


3 


3 


M(Dg K + ) efficiency dependence 


7 


9 


Total 


22 


22 



5 RESULTS 

Figure |21 shows the niES distributions for the reconstructed candidates B~ — * DfK~n~ and B~ — > 
D* + K~k~ . For each mode, we perform an extended unbinned ML fit to the mES distributions 
using the candidates from all decay modes combined. We fit the w-es distributions with the sum 
of the function /(tties) characterizing the combinatorial background and a Gaussian function to 
describe the signal. The mean and width of the Gaussian function, the threshold shape parameter 
£, and the numbers of signal (n S i g ) and background (ribkg) events are free parameters of the fit. 
The likelihood function is given by: 

i=l 

where P^ 19 and P\ kg are the probability density functions for the signal and background, N is the 
total number of events in the fit and i is the index over all events in the fit. 

The fit yields 393 ±25 events in the B~ — > DfK~ir~ mode. Taking into acccount the estimated 
peaking background contribution, we obtain 370 ± 26 signal events for B~ — > DfK~ir~ . The 
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Figure 3: tties spectra for the B~ — * DfK~n~ (left) and B~ — > D* + K~ir~ (right) using the 
data. Solid curves show the fit results, as explained in the text. Dashed lines in the signal regions 
correspond to the background components of the fit. 

number of B~ — > D* + K~tt~ signal events from the fit is 164 ± 17 (no peaking contribution is 
subtracted in this mode as it was estimated to be consistent with 0). We also fit signal yields in 
each of the Df modes, fixing the width and central values of the signal Gaussians to that of the 
combined fit, letting the background level and shape float. The ratio of the signal yields between 
submodes is consistent with the expectations from MC. 

The total signal yield in each B~ decay mode is calculated as a sum over modes (i = (f)ir + , 
K*°K + , KgK + ) and is related to the B~ branching fraction B using the following expression: 

n sig = B ■ N BB ■ } \Bj ■ ei, 

i 

where N BB is the number of produced BB pairs, Bi is the product of the intermediate branching 
ratios and ei is the reconstruction efficiency (from Tabled- As an input to the calculation, we used 
branching fraction numbers from 9 and 12J. The relative systematic uncertainties are converted 
into absolute numbers using the measured central values. The results are: 

B{B~ DfK-ir~) = (1.88 ± 0.13 ± 0.41) • 10" 4 

B(B- -> D* s + K-tt~) = (1.84 ± 0.19 ± 0.40) • 10~ 4 

In summary, two decay modes of charged B mesons are observed for the first time. The 
significance of the observation is 14.2cr for B~ — > DfK~n~ and 9.6a for B~ — > D* + K~ir~ . 
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